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ABSTRACT 



^ |. Context. A detailed radiative transfer code has been previously used to model circumstellar ortho-HaO line emission towards six M- 

I ' type asymptotic giant branch stars using Infrared Space Observatory Long Wavelength Spectrometer data. Collisional and radiative 

O ' excitation, including the V2 = 1 state, was considered. 

Aims. Spectrally resolved circumstellar H20(lio-loi) lines have been obtained towards three M-type AGB stars using the Odin satel- 
C/3 . lite. This provides additional strong constrains on the properties of circumstellar HoO, in particular on the chemistry in the stellar 

^ , ' atmosphere, and the photodissociation in the outer envelope. 

Methods. Infrared Space Observatory and Odin satellite H2O line data are used as constraints for radiative transfer models. Special 
consideration is taken to the spectrally resolved Odin line profiles, and the elfect of excitation to the first excited vibrational states of 
^ , the stretching modes (vi = 1 and vj = 1 ) on the derived abundances is estimated. A non-local, radiative transfer code based on the 

QQ ' accelerated lambda iteration formalism is used. A statistical analysis is performed to determine the best-fit models. 

' Results. The H2O abundance estimates are in agreement with previous estimates. The inclusion of the Odin data sets stronger con- 

straints on the size of the H2 O envelope. The H20( 1 lo-l 01 ) line profiles require a significant reduction in expansion velocity compared 
to the terminal gas expansion velocity determined in models of CO radio line emission, indicating that the H2O emission lines probe 
a region where the wind is still being accelerated. Including the V3 = 1 state significantly lowers the estimated abundances for the 
' low-mass-loss-rate objects. This shows the importance of detailed modelling, in particular the details of the infrared spectrum in the 

range 3 to 6/jm, to estimate accurate circumstellar H2O abundances. 
00 , Conclusions. Spectrally resolved circumstellar H2O emission lines are important probes of the physics and chemistry in the inner 

' regions of circumstellar envelopes around asymptotic giant branch stars. Predictions for H2O emission lines in the spectral range of 

the upcoming Herschel/HIFI mission indicate that these observations will be very important in this context. 
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1. Introduction C/O ratios that are < 1, resulting in a dominance of oxygen- 
bearing molecules (Olofsson 120031 and references therein), as 

The asymptotic giant branch (AGB) marks the final stage in the ^^^^ ^j^^ ^^^^^^ ^^^^^ CO. Carbon stars (C-stars), on the 

hfe of low- to intermediate-mass stars 0.8 - 10 M^). During ^^j^^^ j^^^^^^ j^^^^ ^^^-^^ ^ ^ therefore produce mainly 

this stage the stars lose a significant fraction of their mass in carbon-bearing molecules. However, molecules typical of an M- 

an intense wind, forming a circumstellar envelope (CSE). After environment have been found in reasonably high abun- 

the AGB, what remains is a C/O core, which illuminates the ^^^^^^ carbon-rich AGB stars and vice versa (e.g. Olofsson 

expandingCSE for a short while, forming a planetary nebula Qggjj Agundez & Cernicharo |2006l Schoier et al. i2006t . 

(Habing [1996)1. indicating a departure from thermodynamic al chemical equilib- 

The mass loss dominates the evolution of a star on the AGB, num 
and the physical conditions and the elemental abundances dictate 

the environment for the chemistry. It is therefore essential to try Observations of water vapour in AGB stars are important, as 

to understand the conditions in the inner CSE of an AGB star, in H2O is one of the major (and possib ly the dominant) coolant in 

order to get a handle on the mass-loss mechanism (in particular the CSEs (Truong-Bach et al. |1999l l. with a large amount of far- 

for M-type stars; Woitke l2006l Hofner & Andersen l2007l l. and infrared (FIR) transitions. H2O is expected to be the most abun- 

the molecule and dust formation (e.g. Cherchneff' 120061 1. dant trace molecule next to CO in M-type AGB stars, and as a 

The chemistry in the wind is set by the relative abundances source of circumstellar O and OH it is important for the chem- 

of carbon and oxygen. 'Oxygen-rich' (M-type) AGB stars have istry in the envelope. Indeed, high amounts of H2O have been 

found in M-type AGB stars (e.g. Justtanont et al. 20051 Maercker 

Send offprint requests to: M. Maercker at al. 120081 hereafter Paper I). Cherchneff ( |2006i l shows that a 
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departure from thermodynamical equilibrium due to shocks in 
the inner winds of M-type and carbon-rich AGB stars results 
in higher abundances of H2O in the inner parts of the CSEs. 
Previously, high amounts of H2O have also been found in one 
C-type AGB st ar (IR C+10216/CW Leo; Melnick et al. 120011 
Hasegawa et al. 120061 ). indicating the need for a source of H2O 
other than thermodynamical equilibrium chemistry. In the case 
of M-stars, the H2O emission is expected to come from the 
warm, high-density innermost layers of the CSE, as supported 
by detailed radiative transfer modelling (e.g. Truong-Bach et 
al. 119991 Maercker et al. 12008b . Hence, it is an important probe 
of the velocity structure in the region where the wind accelerates 
(Bains et al. 120031 ). 

In Paper I we determined the circumstellar ortho-H20 abun- 
dance in six M-type AGB stars, using the accelerated lambda 
iteration (ALI) technique for the radiative transfer, based on ISO 
LWS spectra between 43 and 197 fim. The ISO observations do 
not provide spectrally resolved line profiles. However, resolving 
the line profiles is essential, not only for the determination of 
accurate abundances, but particularly to gain information on the 
velocity field in the H2O line-emitting region. In this paper we 
present the results of modelling the line profile of the ground- 
state ortho-H20 line at 557 GHz, Iio-loi^ observed with Odin 
(Nordh et al. l2003l l. towards R Cas, R Dor and W Hya. The mod- 
els are based on the results from Paper I, but we now take special 
consideration to the observed line profile, and examine how this 
helps to further constrain the ortho-H20 abundance distribution 
and the velocity structure in the inner CSE. 

Modelling of H2O is not without difficulties. The high optical 
depths and subthermal excitation of H2O complicate the radia- 
tive transfer considerably. Collisional excitation of the ground 
state can be important in particular for high mass-loss rate ob- 
jects (Paper I), but is rather uncertain due to uncertainties in the 
calculated collisional cross sections (Phillips et al. 119961 ). The 
effect of cooling due to H2O line emission in the inner CSE 
adds additional complications. Finally, pumping through the vi- 
brationally excited states due to IR photons at 2.6 /im and 6 fim 
from the central star and circumstellar dust needs to be taken into 
account. 

In Paper I we showed the ability of our ALI code to han- 
dle the optical depths and subthermal excitation of the H2O 
lines. We further estimated the effects of collisions and exci- 
tation through the first excited vibrational state of the bending 
mode (v2 = 1, at 6.3 yum) on the models, and made a crude esti- 
mate of line cooling by H2O. The first excited vibrational states 
of the symmetric and asymmetric stretching modes (vi - 1 and 
V3 = 1, corresponding to excitation at 2.7 //m and 2.6 fim, re- 
spectively) may be particularly important for low mass-loss-rate 
objects. The high absorption in the stellar photosphere due to 
H2O at these wavelengths makes the proper treatment of radi- 
ation at these wavelenghts especially complicated. We present 
here a first order estimate of how this affects the excitation. 

We give a brief introduction to the Odin and ISO observa- 
tions and basic information on the three objects in Section [21 
and an overview of the modelling procedure in Section [3] In 
Sect. |4| we present the results of the modelling based on both 
Odin and ISO data, and discuss their eff'ects on the predictions 
for the upcoming Herschel/HIFI mission. In Sect. |5] we discuss 
the results and to what degree the resolved line profiles can put 
constraints on the size of the H2O line-emitting region and the 
velocity structure. Our conclusion are presented in section|6l 



Table 1. Summary of Odin observations. 



Source 


Orbits 


rim 

[hours] 


/obs 

[Kkms-'J 


Tmb 

[K] 


Epoch 


RCas 


123 


60.8 


0.5±0.1 


0.04±0.01 


2006.9 


RDor 


118 


36.7 


0.6±0.1 


0.13±0.02 


2006.3 












2006.4 


WHya 


185 


53.4 


0.4±0.1 


0.06±0.02 


2002.9 












2003.6 



2. Observations 

2.1. Odin data 

The H20(I iQ-loi) lines (at 557 GHz) observed by Odin towards 
R Cas, R Dor and W Hya are presented in Fig. [Tl The W Hya 
data was first presented in Justtanont et al. ( 120051 ). Table[T]gives 
a summary of the Odin observations, including the number of 
orbits and the epoch of observations. The integrated intensities 
(/obs) are given in Jmb scale, i.e. they are corrected for the main- 
beam efficiency of 90%. The absolute calibration uncertainty 
is 10%. The observations were done using the position-switch 
mode and the autocorrelation spectrometer (Frisk et al. 120031 
Olberg et al. 120031 ). The internal hot load was used for flux cal- 
ibration. The data was rebinned after a sinusoidal baseline was 
fitted and subtracted. The frequency resolution is 1 MHz, which 
corresponds to 0.54 km s ' at 557 GHz. 

2.2. ISO observations 

In Paper I the models were fitted to ISO LWS spectra in the 
wavelength range of 43 - 197//m. The integrated intensities of 
the ortho-H20 emission lines visible in the spectra were mea- 
sured by fitting a Gaussian to the baseline-subtracted data after 
correcting the spectra for cosmic rays. The dominant error in the 
measured intensities is due to the subtraction of a baseline, and 
it is of the order of ~ 30%. For a more detailed description of 
the ISO observations, the data treatment, and the integrated in- 
tensities for our three objects, see Paper I. 

2.3. Observed sources 

Our sources were studied already in Paper I, and we adopt the 
same stellar and circumstellar parameters in this paper. A brief 
summary of the source characteristics is given below. 

2.3.1. R Cassiopeiae 

R Cas is an M-type Mira variable with a period of 431 days 
(Kholopov et al. 119991 ). A period-luminosity relationship based 
on bolometric magnitudes (Feast et al. 119891 ) gives a luminosity 
of 8725 Lq. From this and the modelling of the spectral energy 
distribution (SED), the stellar distance is estimated to be 172 pc. 
The mass-loss rate, determined through modelling of CO radio 
rotational lines in Paper I, is 9 x lO^^Moyr"'. The gas expan- 
sion velocity of the CSE, derived in the CO emission line mod- 
els, is 10.5 km s '. Based on the ISO observations, the ortho- 
H2O fractional abundance is estimated to be 3.5 x lO"'^ (relative 
to H2; hereafter referred to as the abundance). See Sect. 13. II and 
Paper I for details on the modelling procedures. 
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Fig. 1. Odin observations of the ground-state transition of ortho-H20, (1 10 - loi), at 557 GHz towards R Cas, R Dor, and W Hya. 
The stellar systemic velocity, vlsr, is indicated by the dot-dashed hne, and the solid lines indicate VLSR±Ve (ve is the terminal gas 
expansion velocity determined in the CO models). 



2.3.2. R Doradus 

R Dor is a semi-regular M-type AGB star with a period of 338 
days (Kholopov et al. ll999l l. The revised Hipparcos distance is 
59 pc (Knapp et al. 120031 ). The CO emission-line modelling in 
Paper I gives a mass-loss rate of 2 x 10"^ Mq yr"', and a gas ex- 
pansion velocity of 6.0 km s ' . The ortho-HaO abundance, based 
on the modelUng of the ISO data (Paper I), is 3 x 10 

2.3.3. W Hydrae 

W Hya is a semi-regular M-type AGB star with a period of 382 
days (Kholopov et al. '1999"), located at a distance of 78 pc (re- 
vised Hipparcos distance, Knapp et al. 120031 ). The mass-loss rate 
is 1 X 10"^ Mo yr"' (see Paper I), and the terminal gas expansion 
velocity is 7.2 km s '. The ortho-H20 abundance determined in 
Paper 1 is 15 x 10 Previous estimates of the H2O abundance of 
this star, using various modelling techniques, agree with our esti- 
mate (e^Barlow et al. |1996t Zubko & Elitzur l2000l Justtanont 
et al. l20n5l l. 



through modelling of CO emission lines using the Monte Carlo 
method [presented in detail in Schoier & Olofsson 120011 and 
Olofsson et al. 120021 . Only the integrated intensities were fit- 
ted. In the three cases discussed here, there is a large difference 
between the predicted CO line profiles and the observed ones. 
The models consistently produce double-peaked line profiles, 
while the observations show single-peaked ones. That is, in the 
models the CO line emission is spatially resolved. In the case of 
W Hya, Justtanont et al. ( I2005I I managed to obtain a better fit to 
the CO line profiles by reducing the envelope size and increasing 
the mass-loss rate. However, a similar approach does not work 
for R Cas and R Dor, indicating deviations from our standard 
circumstellar model (although the distance uncertainty also has 
an effect on the line profile). For the moment we do not have 
enough observational constraints to warrant a better mass-loss 
rate determination. The mass-loss-rate uncertainty contributes to 
the uncertainty in the absolute values of the derived abundances. 
However, the sensitivity of the H2O line profiles to the different 
adopted parameters is largely unaffected by this mass-loss-rate 
uncertainty. 



3. The circumstellar emission modelling 

The modelling of the circumstellar line emission was done in the 
same way as in Paper I, using radiative transfer codes based on 
the Monte Carlo and ALI methods to model lines from CO and 
ortho-H20, respectively. Both codes include the radiation con- 
tribution from the circumstellar dust, which was modelled using 
the continuum radiative transfer code Dusty (Ivezic et al. 1999). 
The same 'standard circumstellar model' as in Paper I is adopted, 
with constant mass-loss rate and expansion velocity. The inner 
radius of the CSE is set to the dust condensation radius deter- 
mined in Dusty. 

3. 1 . Dust and CO line emission modeiling 

The ortho-H20 line models in this paper are based on the re- 
sults of the dust and CO line emission modelling in Paper I. 
For model details and modelling procedure, see Paper I and ref- 
erences therein. The radiation from the dust is determined us- 
ing Dusty to fit the SEDs, constructed using 2MASS and IRAS 
fluxes. The mass-loss rate (M), gas expansion velocity (ve), and 
radial kinetic temperature profile of the CSE are determined 



3.2. H2O line emission modelling 

We used the same ALI method as in Paper I in order to model 
the ortho-H20 emission lines (model details and modelling pro- 
cedure are described in Paper I). The molecular data for the 
ground state of ortho-H20 was taken from the Leiden Atomic 
and Molecular DatabastQ (LAMDA) (Schoier et al. 2005). The 
data for the vibrationally excited lines was taken from the 
high-resolution transmission (HITRAN) molecular absorption 
database (Rothman et al. |2005l l. The 45 lowest levels in the 
ground state of ortho-H20 and of the first excited vibrational 
state of the bending mode (v2 = 1) are included, and the eff'ect 
of excitation to the first excited vibrational state of the asym- 
metric stretching mode (V3 - 1) is tested (see Sect. I4.5l l. Since 
the Einstein A-coefficients of the symmetric stretching mode 
(vi = 1 - 0) are an order of magnitude lower than for the asym- 
metric stretching, this state is not likely to affect the models and 
will not be considered here. Collisional rates in the ground state 
are taken from collisions between He and H2O, corrected by a 
factor of 1 .4 to account for collisions with H2 (Green et al. |1993b . 



' available at |http://www.strw.leidenuniv.nI/~moldata| 
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Table 2. Results of the best-fit models from paper I (BFiso) and the best-fit models taking the Odin line profile into account 
(BFodin+iso) for R Cas, R Dor and W Hya. Numbers within brackets give the best-fit values when excitation through the V3 = 1 state 
is included. 



Source BFiso BFo<ji„+iso 





M 






fo 


i5odiii 


<5iso'' 






''e 


fo 


C'Odin 


OISO 






[lO-'Moyr-'] 


[kms-'] 


[10'^ cm] 


[10-»] 


% 


% 




[kms-'] 


[10'^ cm] 


[10-'] 


% 


% 




RCas 


9 


10.5 


3.0 


3.5 


+3 


-5 


0.3 


8.5 


3.0 


3.5 


1.2 


-3.6 


0.3 




















(3.0) 


(3.0) 


(6.8) 


(-1.6) 


(0.3) 


RDor 


2 


6.0 


1.1 


3.0 


+46 


-1 


0.6 


4.0 


1.1 


1.5 


4.5 


13.0 


0.6 




















(1.3) 


(0.6) 


(-3.7) 


(-5.8) 


(0.6) 


WHya 


1 


7.2 


2.2 


15.0 


+67 


-1 


0.5 


5.4 


1.5 


15.0 


-0.3 


-7.7 


0.5 




















(2.2) 


(2.0) 


(-3.0) 


(-7.8) 


(0.4) 



" I'e is determined by the CO line modelling in the BFjso models, and by the Odin line profile in the BFodin+iso models. 
^ The dilference 6 between model and observed intensities is an average difference of all lines for the ISO lines. 



Collisions within the V2 - I state are taken to be the same as in 
the ground state, while collisions between the two vibrational 
states are scaled by a factor of 0.01. This scaling factor lies 
between the ones used by Menten et al. ( |2006| ) and Josselin et 
al. ([2007 ) (0.02 and 0.005, respectively). The dust optical depth 
and temperature profile as determined by Dusty are used. The 
H2O abundance distribution is assumed to be a Gaussian (cen- 
tred on the star) defined by the e-folding radius (where the 
abundance has decreased to 37% of the initial value, fo): 









- 


ir 


n 


V 







The mass-loss rate, gas expansion velocity, and kinetic tem- 
perature profile are included from the CO line emission mod- 
elling. An estimate of is given by the H2O photodissociation 
calculations of Netzer & Knapp (11987b . We present 'best-fit' 
models where the ortho-H20 abundance (fo) and the size of the 
abundance distribution (in terms of the e-folding radius re) were 
treated as free parameters, making it possible to set constraints 
on the photodissociation in the outer envelope (by constraining 
the size of the envelope), and the chemistry in the innermost en- 
velope (by constraining the initial value of ortho-H20). The best 
model fits are determined by fitting the model line fluxes to ob- 
served fluxes using the x' statistics. Paper I gives a detailed de- 
scription of the numerical implementation of the ALI technique, 
as well as of the dependence of the model results on adopted 
numerical and physical parameters of the CSEs. 

Basic results from the dust, CO line, and ortho-H20 line 
emission modelling are given in Table [2l including the reduced 

values of the models, x^^ij The error in the determined abun- 
dances is estimated to be 50% within the circumstellar model, 
and within a factor of ~ 5 in absolute uncertainty. 



4. Results 

While Paper I concerns the determination of the ortho-H20 
abundance based on the spectrally unresolved ISO LWS lines, 
we here model the line profile of the ortho-H20 557 GHz line 
observed with the Odin satellite in order to further constrain the 
models and derive new results on e.g. abundance, molecular en- 
velope size, and gas expansion velocity (in the region where 
the 557 GHz line is emitted). The effect of these new results 
on the line profiles for lines within the range of the upcoming 
Herschel/HIFI mission is also examined. 



4.1. The 557 GHz Odin line predicted by tine ISO models 

Figure [2] shows the H2O line at 557 GHz observed towards 
R Cas, R Dor and W Hya, overlaid by the model lines based 
on the BFiso models from Paper I (dashed line). The differences 
between model and observed fluxes for the Odin and ISO lines 
(^odin and ^iso, respectively) are given in % in Table|2] Most no- 
tably, the model lines are broader than the observed lines. This 
indicates that the 557 GHz line probes a region with a lower ex- 
pansion velocity than the CO lines. In addition, the integrated 
flux of the Odin line is somewhat overestimated in all cases. Of 
importance here is also the fact that, for all three objects, the 557 
GHz line is optically thick throughout the H2O envelope, hence 
the emitting region of this line is limited by the envelope size, 
whereas the ISO line regions are excitation limited. This makes 
the 557 GHz line particularly sensitive to the outer radius of the 
H2O envelope. 

4.2. Line shape dependence on CSE parameters 

Test models were run to determine the general effect of indi- 
vidual parameters on the H20(lio-loi) line profile. The main 
parameters which may affect the line profile are the abundance 
fo, the size of the envelope r^, the shape of the abundance distri- 
bution, the turbulent velocity width Vturb, and the expansion ve- 
locity Ve (which, at a constant mass-loss rate, effectively changes 
the density profile). As an example. Fig. |3] shows the effects of 
these different parameters on the line profile of R Dor. The dif- 
ference in the line fluxes for the Odin and ISO lines compared to 
the BFiso model fluxes is given in Table[3] Other parameters (the 
mass-loss rate, kinetic temperature profile, collisional excitation 
efficiency, the dust and stellar radiation fields, and the excitation 
to the V2 state) affect the integrated intensity of the 557 GHz line 
to the same degree as the ISO lines (tested in Paper I), while its 
line profile does not change significantly. 

In terms of the total flux, Table[3]shows that the 557 GHz line 
is sensitive to the ortho-H20 abundance and outer radius of the 
H2O envelope (in particular compared to the sensitivity of the 
ISO lines to the outer radius). Whereas the integrated line flux is 
less sensitive to the remaining parameters, the shape of the line 
changes noticeably. In particular, the reduction in expansion ve- 
locity results in a narrower line profile and a significantly higher 
peak intensity (+23% in the case of R Dor compared to the best- 
fit ISO model, see Fig.[3j. Changing the abundance distribution 
(/-distr) to a step function (i.e. a constant abundance fo, which 
suddenly drops to zero at r^) results in less self-absorption in the 
blue wing of the line, and it leads to a decrease in the peak inten- 
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Fig. 2. The ortho-H20 (lio - loi) line at 557 GHz observed with Odin towards R Cas, R Dor, and W Hya, overlaid by the BFiso 
model from Paper I (dashed line), and the BFodin+iso model (solid line). The stellar systemic velocity, vlsr, is indicated by the 
dot-dashed Une, and the solid lines indicate VLSR±Ve (ve is the gas terminal expansion velocity determined in the CO model). 




-10 10 
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Fig. 3. The H20(l lo-loi) line profile dependence on various pa- 
rameters. One parameter is varied at a time (the varied parameter 
is indicated in the left part of the figure). ISO refers to the best-fit 
model based on ISO data only from Paper I. The label const prof 
refers to the model with a step function abundance distribution. 
The parameters are varied according to Table [3] (only the results 
of a reduction in abundance and radius is shown). The stellar 
systemic velocity vlsr is indicated by the dot-dashed line, and 
the solid lines indicate VLSR±Ve (ve is the terminal gas expansion 
velocity determined in the CO models). 



sity (-11% for R Dor). The reduced self-absorption in models 
with a step function abundance distribution is due to the sudden 
cut-ofF of the ortho-H20 abundance at r^, removing H2O present 
at further distances responsible for the increased self absorption 
in models using a Gaussian distribution. Increasing the turbulent 
velocity width increases the amount of self-absorption. As can 
be seen in Fig.|3] the most significant change in the shape of the 
line profile occurs when changing the abundance distribution. 

Note that the models predict a significant amount of emission 
in the line wings at velocities outside the range VLSR±i'e- Due 
to the relatively small size of the H2O CSEs, lines with optical 
depths around or larger than unity in the outer parts of the CSEs 
are expected to have a significant contribution of the local line 
profile to the line wings. This means that the observed 557 GHz 
lines are expected to be wider than 2ve. 



Table 3. The change in integrated flux compared to the best-fit 
models using ISO data only. Only one parameter is varied at a 
time. 



Param. 


change 


RCas 


RDor 


WHya 






<5oiiin 


5iso 


"Sodin 


i5iso 


<5odin 


(Jiso 


fo 


+ 100% 


+30% 


+21% 


+ 19% 


+22% 


-2% 


+ 15%, 




-50% 


-20% 


-21% 


-6%, 


-18% 


-11% 


-12% 




+ 100% 


+50%) 


+ 11% 


+30% 


+ 11% 


+34%, 


0% 




-50% 


-35%, 


-14%, 


-32% 


-18% 


-34% 


-18% 


/-distr. 


step func. 


+7% 


+7% 


+ 1% 


+9% 


-6% 


+9% 


V'e 


-20%. 


+6% 


+9% 


^% 


+8% 


-9% 


+3% 


^'turb 


+50% 


+2% 


+ 1% 


+7% 


+9% 


+7% 


+4%, 


yj = 1 


incl. 


+ 10% 


+ 17% 


+2% 


+35% 


+65% 


+24%. 



4.3. Od/n and /SO dafa combined 

In Paper I, the ortho-H20 models were fitted using only the ISO 
LWS spectra, giving rather uncertain constraints on the ortho- 
H2O abundance and, in particular, the size of the H2O abundance 
distribution. The Odin line at 557 GHz is well excited through- 
out the H2O envelope, and it is therefore likely to set a tighter 
constraint on especially the spatial extent of the H2O envelope. 
The eff'ect of the 557 GHz line is tested by fitting models in- 
cluding both the ISO and Odin integrated fluxes. The models are 
varied in ortho-H20 abundance and outer radius r^, and a new 
analysis is done. 

The gas expansion velocity is adjusted to fit the line widths 
of the H20(lio-loi) lines. This requires a reduction in Vg, com- 
pared to that obtained from the CO models, by 19%, 33%, and 
25% for R Cas, R Dor and W Hya, respectively. Although the 
data has only moderate signal-to-noise ratios, a lower expansion 
velocity not only produces a significantly better fit to the line 
shape, but also a better fit to the integrated flux of the Odin line 
and lower values in all three cases. Furthermore, observations 
of the para-H20 line at 183 GHz in maser emission (Gonzalez- 
Alfonso et al. |I998l l indicate expansion velocities for W Hya and 
R Cas of X! 5 km s"' and w 7 km s ', respectively. Thus, we are 
confident that the lower gas expansion velocities derived from 
the Odin lines are real. 

Figure |4] shows the resulting j^'^ maps. Compared to the cor- 
responding maps presented in Paper I (Fig. 3 in Paper I), it is ob- 
vious that the inclusion of the 557 GHz line sets much stronger 
constraints on the envelope size than models only including the 
ISO data. For R Cas and W Hya the models with the lowest 
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values also give the best-fit to the 557 GHz line profile. A slight 
reduction of the abundance by 50%, compared to the lowest 
model, gives an even better fit to the line profile for R Dor The 
best-fit models, BFodin+iso, are marked with triangles in Fig. |4] 
and their results are given in Table |2] For R Dor the model with 
the lowest value is marked by a dot in Fig.|4l 

The resulting 557 GHz line profiles are shown in Fig. |2]as 
the solid lines. The difference between the model ISO and Odin 
fluxes and the observed fluxes is given in % in Table |2] (for the 
ISO lines the average difference is given). In all cases the differ- 
ence is well below the observational error. For R Cas and R Dor 
the radius is the same as in the BFiso models - the inclusion 
of the 557 GHz line only leads to a tighter constraint on the en- 
velope size. For W Hya the inclusion of the 557 GHz line results 
in a new best-fit radius that is 32% smaller, compared to the 
BFiso model. The best-fit abundances do not change for R Cas 
and W Hya, while inclusion of the 557 GHz line results in a 
decrease of the abundance of 50% for R Dor. 



4.4. H/F/ predictions revisited 

Figure|5]shows the difference in the HIFI line profiles for W Hya 
between the BFiso models (dashed lines) and the BFodin+iso 
models taking the 557 GHz line profile into consideration (solid 
lines). The decrease of the outer radius causes the HIFI lines to 
decrease somewhat in integrated intensity, whereas the reduced 
expansion velocity leads to an increased line intensity. These ef- 
fects generally cancel each other in the optically thick lines. For 
optically thin lines (e.g. the two masering lines 532 - 441 and 
634 - 541, and the lines 854 - 76i, 854 - 927, and 845 - 752) the 
reduction in expansion velocity results in an increase of the in- 
tegrated and peak intensities of approximately a factor of two. 
Thus, the prediction in Paper I that these lines are observable 
(some of them easily) with Herschel/HIFI still holds. Note that 
the HIFI lines come from different parts of the CSE and differ- 
ent expansion velocities may apply to them. This will of course 
affect the line profile and results presented in Fig. |5] should be 
regarded with some caution. 

In terms of the line profile dependence on CSE parameters, 
the HIFI lines react similar to the Odin line in Fig.[3j i.e. reduced 
self-absorbtion in step function abundance distribution models, 
narrower lines with reduced expansion velocity, and decreased 



peak and integrated intensities with a decrease in abundance 
and/or outer radius. 



4.5. Effects of excitation tinrougin tine V3 - 1 state 

In order to estimate the effect of excitation through the V3 = 1 
state, we include the 43 lowest energy levels of this state in our 
model and repeat the x~ analysis. The new best-fit results are 
presented in brackets in Table |2] The change in integrated line 
fluxes when including the V3 state are presented in Table [3] The 
average change in integrated flux for the ISO lines is dominated 
by the 743 - 643 transition at 1 16.78 //m for R Dor and W Hya 
(the line changes by 170% and 150%, respectively). Excluding 
this line results in an average change in the ISO lines of 25% and 
15% for R Dor and W Hya, respectively. 

For all three objects the change in the H20(l lo-loi) line pro- 
file is very small when including the V3 = 1 state. A decrease in 
optical depth leads to less self absorption on the blue side and 
less broadening on the red side. 

The Odin line in the W Hya models is, perhaps somewhat 
surprisingly, affected more strongly than in the other two objects. 
This is due to the combination of a low mass-loss rate and the 
very high initial abundance in the best-fit model when only the 
V2 state is included. [The change in the integrated flux using an 
abundance one order of magnitude lower is +4% for the Odin 
line and ^ 40% on average for the ISO lines, more in line with 
the other objects. Also, increasing the mass-loss rate by a factor 
of two moderates the effect of the V3 = 1 excitation.] 

Inspecting the effect of excitation through the V3 state on 
the HIFI lines we find that the integrated flux increases on av- 
erage by 25%, 54%, and 58% for R Cas, R Dor, and W Hya, 
respectively. However, a few lines change by considerably more 
(w 100 - 200%), while others are much less sensitive (with 
changes in integrated flux of only a few percent). 

It is clear that a detailed sampling of the radiation from the 
stellar photosphere at 2.6 jjm is required to properly include the 
first excited vibrational state of the asymmetric stretching mode 
(v3 = 1) in the models. Our models approximate the stellar ra- 
diation by a simple black-body, likely overestimating the flux at 
2.6 /^m (see Sect. 15. 5t . 
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5. Discussion 

5.1. Additional constraints 



of the H2O envelopes determined for R Cas and W Hya in this 
paper. 



The Odin line profiles can be fitted simultaneously with the ISO 
lines without requiring a change in the simple approximations 
made in the 'standard circumstellar model'. However, since we 
are only fitting one resolved line, the assumptions in the cir- 
cumstellar model (e.g. constant Ve, constant M) may be correct 
only within the region from which the line is coming. When 
Herschel/HIFI data becomes available and the line profiles of 
a series of lines have to be fitted, the situation may change. As 
the lines probe different regions in the CSE, fitting different lines 
may require a gradient in the velocity distribution, i.e. an acceler- 
ation zone, and a different abundance distribution (see Sect. 15.41 ). 

5.2. The size of the H2O envelope 

For all three objects, the Odin line at 557 GHz puts an addi- 
tional constraint on the H2O envelope size. Figure |6] shows the 
tangential optical depth of the Odin line for R Cas, R Dor and 
W Hya. The line is excited throughout the envelope, until the 
ortho-H20 abundance drops at the e-folding radius re. In Paper 
I, the BFiso model radius was found to be within 1 cr of the theo- 
retical radius given by Netzer & Knapp ( 119871 1. The latter is ob- 
tained from models of H2O photodissociation and it is defined 
as the distance where the abundance of OH is highest. Although 
the 557 GHz data requires a reduction in the outer radius, the 
new best-fit radii based on the combined Odin and ISO data are 
still within the Icr contours presented in Paper I, though with 
much tighter constraints, confirming the photodissociation mod- 
els of Netzer & Knapp ( |1987t . Observations of the 1665 and 
1667 MHz OH masers in the CSE around R Cas (Chapman et 
al. 1994) and W Hya (Szymczak et al. T99S) show OH masering 
shells at (1.3 - 2.6) x lO'^'cm and (0.8 - 1.3) x lO''^ cm from the 
central star, respectively. This is consistent with the smaller radii 



5.3. The abundance of H2O 

The derived total (para + ortho) H2O abundances (assuming 
an ortho-to-para ratio of 1) in Paper I are ~ 6 xlO The ther- 
mal equilibrium abundance of H2O is approximately 1-3 x lO-** 
(Mamon et al. 119871 Cherchneff 120061 ). while the upper limit 
given by the cosmic abundances of carbon and oxygen is no 
more than 1 x 10"^, assuming that all carbon goes into CO and the 
remaining oxygen into H2O (Anders & Grevesse |1989l ). Models 
of departures from thermodynamical equilibrium chemistry due 
to changes in density and temperature as a result of shocks in the 
inner parts of the winds predict total (para + ortho) H2O abun- 
dances on the order of ^ 4 x 10 "* (Cherchneff l2006b . This likely 
explains at least part of the observed high abundances of H2O. 
Other processes which may explain an increased abundance of 
H2O are the evaporation of icy bodies (Justtanont et al. 120051 ) or 
grain surface reactions (Willacy 12004b . However, the inclusion 
of the Odin line and the V3 = 1 mode in the excitation analy- 
sis result in lower abundances. Considering the current absolute 
uncertainties (at least on the order of a factor of 5), the results 
cannot conclusively confirm or rule out the various competing 
scenarios, from thermodynamical equilibrium chemistry to the 
evaporation of icy bodies. 

5.4. The acceleration zone 

In dust-driven winds, the dust accelerates and initiates the mass 
loss relatively close to the star. Depending on the size of the 
H2O line-emitting region, observations of water vapour lines 
may probe the inner structure of the CSE and give informa- 
tion on the velocity field. A number of observations indicate 
that the H2O line-emitting region coincides with the accelera- 
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tion zone in M-type AGB stars. Imaging in the 22 GHz H2O 
maser line, indicates that the terminal velocity (i.e. the final gas 
expansion velocity) is reached at relatively large distances from 
the central stars (Bains et al. 120031 ). For R Cas, the analysis 
of H2O masers by Colomer et al. ( I2000I I indicates an expan- 
sion velocity of only 3.0 km s ' in the inner regions of the en- 
velope (^ 3.4 X lO'^^cm; the gas terminal expansion velocity 
estimated from CO data is 10.5 km s"'). Similarity, models fitted 
to OH, H2O and SiO masers in W Hya indicate that the terminal 
expansion velocity (estimated from CO data to be 7.2 km s"') 
is reached at distances > 2 x 10'^ cm (Szymczak et al. 119981 ). 
The H2O maser is located in a ring-like structure with a ra- 
dius of 1.8 X 10'"* cm expanding with a velocity of 3.6kms"' 
(Reid & Menten [T990] l, and the OH masers lie at 8.2 x 10"* cm 
and 1.3 x 10'^ cm with expansion velocities of 4.2kms ' and 
5.2kms"\ respectively (Szymczak et al. 1998). The SiO maser 
has its origin at a distance of 5.8 x 10'^ cm (Miyoshi et al. 119941 ), 
and the lower limit of the expansion velocity is estimated to be 
about 2kms ' (Szymczak et al. |1998t . Likewise, in order to 
model high angular resolution observations of SiO in the CSE 
around R Dor, Schoier et al. ( |2004| i determined a best-fit model 
to the data using a two component density distribution (with a 
higher SiO abundance within 1.0 x lO'^'cm) and a lower ex- 
pansion velocity in the inner regions of the CSE of 4.9 km s ' 
(the gas terminal expansion velocity estimated from CO data is 
6.0kms-i). 

The fits to the Odin line profile indicate that the H2O line- 
emitting region covers an area in which the expansion velocity 
has not yet reached its terminal value. As can be seen in Fig.|6l 
the dominant part of the H2O 557 GHz line emission comes from 
a region at a radius of about 10'^ cm for all three sources. For 
R Cas, R Dor, and W Hya, good fits to the 557 GHz line profiles 
require expansion velocities of only 8.5, 4.0, and 5.4 km s"', re- 
spectively. This means a 19, 33, and 25%, respectively, decrease 
from the terminal gas expansion velocities as measured in the 
CO radio lines. Thus, the 557 GHz line data is consistent with 
an acceleration zone extending beyond 10''' cm for these, rela- 
tively, low-mass-loss-rate objects. 

Fig. I2] shows the estimated positions of the masers and the 
corresponding gas expansion velocities for W Hya. The expan- 
sion velocity required to fit the Odin 557 GHz line profile has 
been placed at a radius of 1 x 10'^ cm (where the brightness dis- 
tribution peaks), and the terminal gas expansion velocity deter- 
mined by the CO emission lines is placed at 1.0 x 10'*'cm (the 
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estimated CO photodissociation radius, Olofsson et al. 12002b . 
The measurements show that the wind has not reached its ter- 
minal value in the region where the 557 GHz line originates. 
Figure [8] shows the tangential optical depth distributions of four 
of the HIFI lines shown in Fig. |5]for W Hya. Two of the lines 
probe a region similar to that of the 557 GHz line, whereas the 
other two appear to be emitted very close to the star. The HIFI 
observations of H2O will likely probe different regions of the 
CSE, and in particular may probe regions with different expan- 
sion velocities. Fitting the individual line profiles of the HIFI 
observations should therefore give important information on the 
velocity structure of the CSE. 

5.5. The importance of detailed models 

Inclusion of excitation to the first excited vibrational state of the 
asymmetric stretching mode (V3 - 1) changes the integrated 
fluxes of the modeled lines significantly, in particular for the 
low-mass-loss-rate objects R Dor and W Hya. On the contrary, 
Gonzalez-Alfonso et al. (12007b . when modelling the lio - loi 
transition observed with the Submillimeter Wave Astronomy 
SatelUte (SWAS) in the carbon star IRC+10216, concluded that 
the lio - loi transition is not significantly affected by pumping 
through the V3 = 1 state, as the pumping of the line occurs mainly 
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through the vt-I state in this object. However, the abundance 
of 0-H2O in their models is low, and IRC+10216 has a high 
mass-loss rate. We find the same trend with mass-loss rate and 
abundance in our models of the lio - loi emission. Models of 
photon dominated regions show that IR pumping affects higher 
lying transitions more strongly in accordance with our results for 
lower abundances (Gonzalez Garcia et al. 120081 ). 

It must be emphasized here that our models approximate 
the radiation from the central star with a blackbody, and there- 
fore likely overestimate the flux at 2.6 fim as photospheric ab- 
sorption by H2O is not taken into account. This is confirmed 
by low-resolution SWS spectra from ISO for R Dor (Ryde & 
Eriksson |2002] | and W Hya (Justtanont et al. 12004b . It is there- 
fore difficult to accurately include the effect of the V3 - 1 state 
without detailed modelling of the photosphere. The H2O abun- 
dance for the low-mass-loss-rate objects will likely lie between 
the abundances determined when including both the V2 = 1 and 
V3 = 1 states and when including only the V2 = 1 state (assum- 
ing a blackbody radiation field). This demonstrates how impor- 
tant detailed models will be in order to determine accurate H2O 
abundances from Herschel/HIFI data. 

6. Conclusions 

We have modelled the ortho-H20(l 10 - loi) line at 557 GHz ob- 
served with the Odin satellite for the M-type AGB stars R Cas, 
R Dor and W Hya. The models are based on previous models of 
ortho-H20 lines observed with ISO (Paper I). The Odin line pro- 
file is resolved in contrast to the ISO lines, and therefore gives 
additional constraints to the circumstellar model. Models that 
fit the integrated intensities of both the ISO and Odin lines re- 
sult in tighter constraints on the size of the H2O envelope than 
models based on ISO observations only. The resulting H2O en- 
velopes are smaller than derived from the ISO data alone. The 
envelope sizes given by models of H2O photodissociation lie 
within the errors of our estimates, confirming the photodisso- 
ciation chemistry in the outer envelope. Additionally, expansion 
velocities lower than those estimated from CO data are required 
to fit the 557 GHz line profile, hence indicating that observations 
of water vapour probe the acceleration zone in the CSE. The es- 
timated sizes of the H2O envelopes and the expansion velocities 
are confirmed by observations of H2O and OH maser emission. 

The H2O abundances from Paper I are confirmed here and, 
taken at face value, cannot be explained by thermodynamical 
equilibrium chemistry. However, the effect of excitation through 
the first excited vibrational state of the asymmetric stretching 
mode (v3 - 1) is tested (assuming a stellar blackbody field), and 
it leads to significant decreases in abundances for the low-mass- 



loss-rate objects. This shows the difficulty in accurately mod- 
elling circumstellar H2O emission lines. Careful and detailed 
modelling of the stellar radiation field is necessary in order to 
determine correct H2O abundances, and hence identify the ori- 
gin of water vapour in the CSEs of AGB stars, and to determine 
important properties of the inner CSEs. This will be particularly 
important for the upcoming Herschel/HIFI mission. 

The HIFI observations will help to set constraints on the H2O 
abundance distribution as well as the kinematics of the CSE. 
HIFI/Herschel will observe H2O emission lines in the range 557 
to 1885 GHz, covering a range of excitation conditions and op- 
tical depths. While the lower lying transitions (like the 557 GHz 
line) will constrain the size of the emitting region (as they are 
sensitive to the outer radius), higher lying transitions will set 
constraints on the velocity structure close to the star. As these 
two parameters are constrained by fitting the line profile for a 
range of emission lines, the abundance estimates made in Paper 
I and here will likely be improved significantly. 
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